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Correlation of the angular dependence of spin-transfer torque and giant 
magnetoresistance in the limit of diffusive transport in spin valves 
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Angular variation of giant magnetoresistance and spin-transfer torque in metallic spin-valve het- 
erostructures is analyzed theoretically in the limit of diffusive transport. It is shown that the 
spin-transfer torque in asymmetric spin valves can vanish in non-collinear magnetic configurations, 
and such a non-standard behavior of the torque is generally associated with a non-monotonic angu- 
lar dependence of the giant magnetoresistance, with a global minimum at a non-collinear magnetic 
configuration. 

PACS numbers: 



Introduction: The macroscopic model of electronic 
transport in magnetic metallic multilayers, developed by 
Valet and Fert , is commonly used for interpretation of 
experimental results on current-perpendicular-to-planc 
(CPP) giant magnetoresistance (GMR). The model in- 
cludes several phenomenological parameters which can 
be extracted from fitting experimental and theoretical re- 
sults on resistance/magnetoresistancc in coUinear (paral- 
lel and antiparallel) magnetic configurations. The model 
is based on the assumption that the spin diffusion length 
is longer than mean free path, and the latter is smaller 
than the layer thickness. Recently, Penn and Stiles Q 
showed that the Valet-Fert model is justified even for 
spin diffusion lengths comparable to the mean- free paths. 
Moreover, the model fits well to experimental results even 
when the mean free paths are comparable to the layer 
thicknesses. 

Validity of the model, however, can be justified only 
when it gives results which are in agreement with exper- 
iment not only in the coUinear configurations (too many 
fitting parameters), but also in the full range of angles be- 
tween magnetic moments of the layers' magnetizations. 
This problem was not addressed yet. In a recent pa- 
per [3| the Valet-Fert model has been extended to non- 
collinear configurations, but only the angular variation of 
the spin-transfer torque (STT) was analyzed there. The 
spin torque plays a crucial role in the phenomenon of 
current-induced magnetic switching (CIMS). It turned 
out that both CPP-GMR and CIMS effects are corre- 
lated. Moreover, there are normal and inverse GMR and 
also normal and inverse CIMS phenomena. Accordingly, 
four different possibilities can be found in real systems 
Indeed, the normal and inverse CPP-GMR and/or 
CIMS have been demonstrated by manipulatiiig the bulk 
and/or interface spin-asymmetry parameters [J]. 

It has been shown in Ref.0] that the STT in asym- 
metric structures can vanish at a noncoUinear configura- 
tion, which has a significant impact on the stability of 
magnetic configuration. As a result, precessional states 
in zero magnetic field were predicted in the Co/Cu/Py 



nanopillars Q and later experimentally confirmed 
One may naturally expect that the nonstandard behav- 
ior of STT may be associated with some anomalous an- 
gular behavior of the CPP-GMR. This problem is ad- 
dressed here, and wc show that a global minimum in 
resistance of asymmetric structures may occur in a non- 
collinear configuration. This non-monotonic behavior of 
the resistance (and consequently also GMR) is generally 
accompanied by the non-standard angular dependence 
of STT. These two features seem to be characteristic of 
asymmetric systems in the diffuse transport regime, pro- 
vided the system's parameters obey some conditions. 

Model of CPP-GMR: Within the diffusive approach 
[II, S 0] I spatial dependence of the average electrochem- 
ical potential in a ferromagnetic (F) layer has the gen- 
eral form; /I = —(3 g -\- Cx -\- G, where the axis x is nor- 
mal to the structure, and g is the spin accumulation, 
g = Aexp{x/lsf) -\- B ex.p{—x/lsf), with Igf being the spin 
diffusion length. Similar formula also holds for normal 
metal (N) layers, but with j3 = 0. All the constants {A, 
B, C, and others) entering the general solution of the 
diffusive equations in different layers can be determined 
from the corresponding boundary conditions [J, S 0] . 

The driving field can be then calculated as E{x) = 
{\/e){dfjL/dx) The presence of N/F interfaces gives 
rise to additional voltage drops due to spin accumulation 
in their vicinity. The total voltage drop can be then writ- 
ten as Ay = X^i^^i! where AVi is the voltage drop in 
the i-ih layer of the spin valve (voltage drops at interface 
resistances will be included to the ferromagnetic layers). 
When the index i corresponds to a ferromagnetic layer, 
Al^i = /SV^^ -\- AV^"^"^. If, however, i corresponds to a 
nonmagnetic layer, AT^i = Al/j^^ -I- pidilo, where pi is 
the bulk resistivity of the normal metal, di is the corre- 
sponding layer thickness, and Iq is the current density. 
The voltage drops due to spin accumulation (in magnetic 
and nonmagnetic layers) read 



AFf = / [E{x)-Eo]dx, 



(1) 
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where the corresponding electric field Eq is taken far 
from the interface. Apart from this. AV^^^ = Iq i?^^' 
for ferromagnetic films, i?-"^' = + {1/Rii)]~^, 

Ria = Rl^ + dp.,„ + Rf^ (for a i), with pi„ being 
the corresponding spin dependent bulk resistivity, and 
i?^ {Rf„) denoting the intcrfacial resistances (per unit 
square) associated with the left (right) interface of the 
i-th (ferromagnetic) film. 

The total resistance of the system (per unit square) 
is i? = Ay/ Jo, while the magnetoresistance, AR{9) = 
R{9) — i?(0), describes a change in the total system re- 
sistance when magnetic configuration varies from a non- 
collinear to parallel one. We note that what one needs 
to calculate are the AV^^ contributions only. It is con- 
venient to define reduced magnetoresistance as 



havior of the reduced GMR is more pronounced in spin 
valves that are more asymmetric, see Figllje). 



R{0) - Rp 
Rap - Rp 



(2) 



Several theoretical approaches have beenproposed to 
describe the angular variation of GMR [8|. The ex- 
plicit form for the reduced magnetoresistance, r{9) = 
sin^(0/2)/[l -|-xcos^(0/2)l has been extracted within the 
magneto-circuit theory Q and diffusive approach [lo| . 
In recent measurements on Py/Cu/Py valves [ll| the 
parameter x has been treated as a fitting parameter, 
and the r(0) dependence was found to describe the ex- 
perimental data relatively well. However, this formula 
breaks down for asymmetric spin valves, where a global 
minimum of the system resistance may appear at a non- 
collinear configuration [llj . 

Results: The minimum in resistance (and also in GMR) 
in a noncoUinear configuration appears only in asymmet- 
ric spin valves. In the following we discuss the angular 
dependence of the STT and GMR in spin valves for pos- 
itive current density, /q > (electrons flow from right 
to left, or in other words charge current flows from the 
layer of thickness di to the layer of thickness d2). Figure 
[Ija) shows electric field profile in the Co(di)Cu(10)Co(8) 
spin valve for di ^ 16 nm sandwiched between semi- 
infinite Cu leads and for both parallel (P) and antiparallel 
(AP) magnetic configurations. The angular dependence 
of the corresponding voltage drops within the Co lay- 
ers is shown in Fig[ljb). The voltage drop within the 
Co(16) layer exhibits a very weak minimum for 9 ~ 7r/3. 
The minimum becomes much more pronounced for larger 
layer thickness, as shown in Fig[T]Jc) for di = 60 nm. 
Since the total voltage drop is a sum of all drops in the 
individual layers, the GMR can exhibit a minimum at a 
non-coUinear configuration when the resistance decrease 
in the thick F layer overcomes the resistance increase in 
the thin F layer. The global minimum arises as a result 
of the spatial depletion of electrical field in the thick F 
layer, which is a consequence of spin accumulation dis- 
continuity at the N/F interface controlled by the mix- 
ing conductances. The reduced GMR, r{9), is shown in 
FiglHd) for both values of di. The non-monotonic be- 
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FIG. 1: (Color online) Transport characteristics of the biased 
Cu/Co(di)/Cu(10)/Co(8)/Cu spin valve, (a) Spatial depen- 
dence of the electric field in the system for di — 16 nm and 
for P and AP configurations. Angular dependence of the volt- 
age drops in Co(di) shown by dashed (blue) line and Co(8) 
by solid (red) line for (b) di — 16 nm and (c) di = 60 nm. 

(d) Angular dependence of the reduced magnetoresistance. 

(e) Reduced magnetoresistance as a function of 6 and di. In 
parts (a), (b) and (c) the current density Jo = 10* A/cm^ was 
assumed. The other parameters are as in Ref. 0]. 

In Figs [2Ka)-(c), the diagrams present the regions of 
layer thicknesses, where the non-monotonic behavior of 
the reduced GMR can be observed (gray regions). For 
the Co/Cu/Co spin valves [Fig[2I^a)] as well for the 
Py/Cu/Py ones [Figl^I^b)], the diagrams are symmetric 
with respect to c?i = ^2, and the non-monotonic angular 
variation of the GMR (global minimum at a non-coUinear 
configuration) can be noticed for spin valves with sig- 
nificantly different layer thicknesses. In Co/Cu/Py spin 
valves, where an additional asymmetry appears due to 
different magnetic materials, a non-monotonic angular 
variation of the GMR can be observed even for compa- 
rable layer thicknesses, see FiglJl^c). This is mainly due 
to strong asymmetry in spin diffusion lengths of Co and 
Py. but difference in the bulk as well as interface spin 
asymmetries of the Co and Py also contributes to the 
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non-monotonic behavior. 
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FIG. 2: (Color online) Diagrams illustrating presence of a 
global magnetoresistance minimum at non-coUinear config- 
urations - gray regions - and angular spin-transfer torque 
dependences. (a) Diagram for the Co(di)Cu(10)Co(ci2), 
(b) Py(di)Cu(10)Py(d2), and (c) Co(di)Cu(10)Py(d2) spin 
valve. The solid (red) and dashed (blue) lines denote critical 
thicknesses where 9r/90|e_,o = for the torque exerted on 
the left F layer of thickness d\ and right layer of thickness d^ , 
respectively. (d,e,f) Angular dependence of the spin-transfer 
torques acting on the left F layer of thickness di - shown by 
solid (red) lines - and right F layer of thickness d2 - shown by 
dashed (blue) lines - for systems corresponding to the dots in 
the left panel (a,b,c). 

Experimental observations on Py/Cu/Py spin valves 
revealed a weak non-monotonic angular variation of the 
GMR effect 
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This has been attributed to the absorp- 
tion of transverse spin accumulation in a non-coUinear 
configuration, which reduces the resistance. Such absorp- 
tion also gives rise to the STT acting on the F layer [l^ , 
which in asymmetric spin valves can exhibit an anoma- 
lous (non-standard) angular dependence 0, @|. In sys- 
tems with a non-standard STT, the transverse compo- 
nent of spin current (accumulation) at the active N/F 
interface vanishes at a certain noncoUinear configuration. 
The presence of a GMR minimum in a noncoUinear con- 
figuration can be thus related to the non-standard STT. 

The STT appears due to absorption of the transverse 
spin current component at the N/F interface [l^ . and 



where the superscripts L and R denote the left and right 
interface, respectively, associated with the F layer. De- 
pendence of the STT can also be expressed explicitly in 
terms of the mixing conductances and spin accumulation 
at the normal-metal side of the N/F interface d, 0]. The 
STT consists generally of both in-plane and out-of-plane 
components. Since the latter component is much smaller 
than the former one (due to small imaginary part of the 
mixing conductances [l^l ) , in the following discussion we 
will consider only the in-plane component. In asymmet- 
ric spin valves, the proper choice of magnetic materials 
and/or layer thicknesses can result in vanishing STT at 
a non-coUinear magnetic configuration Such a non- 
standard STT destabilizes both coUinear configurations 
for positive current and stabilizes both configurations for 
negative current [^, The former case is of particu- 
lar interest as the non-standard torque leads to current- 
induced steady state oscillations in the absence of exter- 
nal magnetic field 0, 0|. In FiglJId) we show the angular 
dependence of STT in the Co(60)Cu(10)Co(8) spin valve 
exerted on the Co(60) (sohd line) and Co(8) (dashed 
line). The STT acting on the Co(8) layer destabilizes 
P and stabilizes AP configuration, whereas the torque 
acting on the Co(60) vanishes at a non-coUinear configu- 
ration and stabilizes both P and AP configurations. The 
torques in Fig^d) correspond to the system indicated 
by the dot in Fig^a). This point is below the critical 
line, given by dr/dOle^o = 0, which identifies the region 
where a non-standard STT acting on the Co{di) layer ap- 
pears. When the layer thicknesses are above the critical 
line, but still in the gray region, see the dot in Figlll^b), 
the torque acting on the particular F layer vanishes only 
for the coUinear configurations, as shown in Figl^fe) for 
the Py(l)Cu(10)Py(6) spin valve, but reduced GMR stiU 
exhibits a global minimum for a non-coUinear configu- 
ration. Since the critical lines are close to the bound- 
ary of the non-monotonic angular GMR behavior (gray 
regions), the non-standard STT is correlated with the 
non-monotonic angular variation of GMR. 

At the critical angle 9c, where the torque r vanishes, 
the transverse component of spin accumulation at the ac- 
tive interface disappears. In a general case, however, the 
angle 9g between the spin moment of the F layer and spin 
accumulation vector at the normal- metal side of the N/F 
interface is nonzero. Angular dependence of the STT 
can be then expressed as a function of 9g In Fig[3]Ja) 
and (b) the in-plane spin accumulation components at 
the normal-metal side (in the spacer layer) at the left 
and right interfaces are shown for the spin valves con- 
sidered in Figs [2l^d)-(f). The components are expressed 
in local coordinate frames, where gz points along mag- 
netization in the left F layer whereas 5^ along magneti- 
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zation in the right F layer. The STT acting on Co(8) in 
the Co(8)Cu(10)Py(8) spin valve exhibits regular behav- 
ior. Spin accumulation in the P configuration is positive 
and has comparable amplitudes in the vicinity of both 
interfaces [see the dots on the dotted lines in Figs [3]Ja) 
and (b)] due to long spin- flip length in Cu (Zsf — 1 A^m) . 
When magnetization of the right layer rotates in the film 
plane, spin accumulation at the left interface roughly fol- 
lows the net-spin of the Py(8) layer, and the angle 9g 
is a monotonic function of 9, see FiglSl^c). Figure [2ljf) 
shows that STT exerted on the Py(8) layer vanishes in 
a non-coUinear configuration {0 = dc), for which the 
component also vanishes. At 6 ~ 9c one finds 9g — 0. 
For the Co(60)Cu(10)Co(8) system, STT acting on the 
Co (60) layer shows a non-standard behavior which is 
qualitatively similar to that for the Py(8) layer in the 
Co(8)Cu(10)Py(8) spin valve. Angular variation of STT 
for the Py(l)Cu(10)Py(6) spin valve, shown in FiglSlJe), 
vanishes regularly in P and AP configurations, and 9g is 
a monotonic function of 9. 
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FIG. 3: (Color online) Parametric plots of the in-plane 
spin accumulation components and angular dependences of 
the angle 6g between the spin accumulation and spin mo- 
ments, (a) spin accumulation at the left (F/N) interface; 
(b) spin accumulation at the right (N/F) calculated at the 
normal metal side (N) in the vicinity of the active interfaces. 
The spin accumulation components are expressed in their lo- 
cal reference frames. The dashed (red) lines correspond to 
Co(60)Cu(10)Co(8); solid (black) lines to Py(l)Cu(10)Py(6); 
and dotted (blue) lines to Co(8)Cu(10)Py(8). The filled dots 
correspond to the parallel configuration of the layer magne- 
tizations. Angular dependence of the 9g at the (c) left F/N 
and (d) right N/F interface. 

Non-collinear configuration of the F layer magnetiza- 
tions leads to discontinuities of the spin accumulation 
at the F/N interfaces [angle 9g in Figs 3(c) and (d)]. 
From this we deduce that if one takes the thickness of 
one of the F layers smaller than the corresponding spin 



diffusion length and thickness for the second F layer is 
larger than the appropriate spin diffusion length, then 
the spin accumulation is predominately determined by 
the later F layer. One finds then non-standard STT and 
non-monotonic GMR angular behavior. We have found 
that this behavior is mostly controlled by the mixing 
conductance of the interface between spacer layer and 
that F layer whose thickness is smaller than the corre- 
sponding spin diffusion length. For instance reducing 
the mixing conductance at the Co(8)/Cu(10) interface 
in the Co(8)Cu(10)Py(8) valve by about 50% lifts the 
non-standard STT and GMR behavior. 

In conclusion, what stems from the above results is 
a need for further experimental investigations, and that 
Co/Cu/Py system is a good candidate to test the theo- 
retical predictions. This also could answer the question 
whether the diffusive approach used to analyze CPP- 
GMR in coUinear configurations is well justified. To ar- 
rive at more convincing conclusions one also should cor- 
relate the results on GMR with those on STT. 
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